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1 Abstract 



The effect of the random magnetic field distribution on the phase diagrams and ground state 
magnetizations of the Ising nanowire is investigated with effective field theory with correlations. 
Trimodal distribution chosen as a random magnetic field distribution. The variation of the phase 
diagrams with that distribution parameters obtained and some interesting results found such as 
reentrant behavior. Also for the trimodal distribution, ground state magnetizations for different 
distribution parameters determined which can be regarded as separate partially ordered phases of 
the system. Keywords: Ising Nanowire; random magnetic field; trimodal distribution 
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Recently there has been growing interest both theoretically and experimentally in the magnetic 
nanomaterials such as nanoparticles, nanorods, nanotubes and nanowires. Nowadays, fabrication of 
these nanomaterials no longer difficult, since development of the experimental techniques permit us 
making materials with a few atoms. For instance acicular magnetic nano elements were fabricated 
[HOE] and magnetization of the nanomaterial has been measured [4] . In general nanoparticle 
systems can be used as sensors [5] and they can be used for making permanent magnets [B], 
beside some medical applications [7]. In particular, magnetic nanowires and nanotubes have many 
applications to a nanotechnology O |9] . Nanowires can be used as an ultrahigh density magnetic 
recording media [101 HH IE] and they have potential applications in biotechnology [TSj [14] , such as 
Ni nanowires can be used for bio seperation [15[ [T6] . 

In a nanometer scale for the magnetic properties of these materials, some unexpected properties 
appear and in a nanotechnology these properties allow us to fabrication materials for various 
purposes, since behaviors of these finite materials different than bulk counterparts. In addition, 
magnetic properties and magnetic phase transition characteristics are highly depends on the size 
^ ' and the dimensionality of the material. 

1 Mean field approximation (MFA), effective field theory (EFT) and Monte Carlo (MC) simula- 
tion are most common used theoretical methods for determining the magnetic properties of these 
materials. For instance, nanoparticles investigated by EFT with correlations y/7], MFA and MC 
[18] . The phase diagrams and the magnetizations of the nanoparticle described by the transverse 
Ising model investigated with using MFA and EFT [T9j [20] , compensation temperature of the nano 
particle [21] and magnetic properties of the nanocube with MC [22] are among these studies. 

Another method namely Variational cumulant expansion (VCE) based on expanding free energy 
in terms of the action up to m*'' order, has been applied to the magnetic superlattices [23] and 
ferromagnetic nanoparticles [24l [25] . The first order expansion within this method gives the results 
^ ; of the MFA. 

, Various nano structures such as FePt and Fe^Oi nanotubes [26] can be modeled by core- 
shell models and these models can be solved also by MFA, EFT and MC. The phase diagrams 
and magnetizations of the transverse Ising nanowire is treated within MFA and EFT [37] [2H], 
the effect of the surface dilution on the magnetic properties of the cylindrical Ising nanowire and 
nanotube has been studied [29l [30], the magnetic propeties of nanotubes of different diameters, 
using armchair or zigzag edges has been investigated with MC [31], initial susceptibility of the 
Ising nanotube and nanowire within the EFT with correlations |32[ 133] and the compensation 
temperature which appears negative core-shell coupling is investigated by EFT for nanowire and 
nanotube [34]. There has been also works deal with hysteresis characteristics of the cylindrical 
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Ising nanowire |351I36| . Beside these, higher spin nanowire or nanotube has been investigated also, 
e.g. spin - 1 nanotube and nanowire |38] . mixed spin - 3/2, 1 core shell structured nanoparticle 
[39] . mixed spin - 1/2, 1 nanotube |40] . 

On the other hand, as far as we know, there is less attention to quenched randomness effects of 
these systems except the site dilution. However, including quenched randomness or disorder effects 
to these systems may induce some beneficial results. For this purpose we investigate the effects of 
the random magnetic field distributions on the phase diagrams of the Ising nanowire within this 
work. As stated in '5^ the phase diagrams of the nanotube and nanowire are similar qualitatively, 
then investigation of the effect of the random magnetic field distribution on the nanowire will give 
hints about the effect of the same distribution on the phase diagrams of the nanotube. 

The Ising model in a quenched random field (RFIM) has been studied over three decades. The 
model which is actually based on the local fields acting on the lattice sites which are taken to be 
random according to a given probability distribution was introduced for the first time by Larkin [41j 
for superconductors and later generalized by Imry and Ma . Beside the equality between the 
diluted antifferromagnets and ferromagnets with random field distribution |431 144j , the importance 
of the random field distribution on these systems comes from the fact that, random distribution of 
the magnetic field changes drastically phase diagrams of the system and then magnetic properties. 

The paper is organized as follows: In Sec. [3] we briefly present the model and formulation. The 
results and discussions are presented in Sec. |4l and finally Sec. [5] contains our conclusions. 

3 Model and Formulation 

We consider a nanowire which has geometry shown in Fig. [T] The Hamiltonian of the nanowire is 



Figure 1: Schematic representation of a cylindrical nanowire (top view). The gray/black circles 
represents of the surface/shell magnetic atoms, respectively 

given by 

H = - Jl ^ SiSj - J2 ^ SmSn - J3 ^ SjSm " ^H^Si - ^i?,„S„ (1) 

<i,j> <m,7i> <i,m> i m 

where Si is the z component of the spin at a lattice site i and it takes the values Si — ±1 for the 
spin-1/2 system, Ji and J2 are the exchange interactions between spins which are located at the 
core and shell, respectively, and J3 is the exchange interaction between the core and shell spins 
which are nearest neighbor to each other. Hi and Hm are the external longitudinal magnetic fields 
at the lattice sites i and m respectively. Magnetic fields are distributed to lattice sites by a given 
probability distribution. The first three summations in Eq. ([1} is over the nearest-neighbor pairs 
of spins, and the other summations are over all the lattice sites. 
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This work deals with the following discrete magnetic field distribution, namely trimodal distri- 
bution 

P {Hi) = p5 {Hi) + [5 {H, -^Ho) + 6 {H, + Ho)] (2) 

which covers a bimodal distribution for p = and reduces to the system with zero magnetic field 
(pure system) for p = 1. According to the distribution given in Eq. ([2]), p percentage of the lattice 
sites are subjected to a magnetic field Hi — 0, while half of the remaining sites are under the 
influence of a field Hi = Hq whereas the field Hi = —Hq acts on the remaining sites. 

We have to separate four representative spins according to their interactions with other spins. 
Their magnetizations {nii^i = 1,2,3,4) can be given by usual EFT equations which are obtained 
by differential operator technique and decoupling approximation (DA) [451 146] , 

mi = [Ai + miB{f [A3 + mzSg] [A3 + m^B^f [Ai + m^Bi] 
1712 = [A-i + miB-i][A2 + 'm2B2f[A2 + m3B2f ^ 
"^3 = [Az + miBz] [A2+ 7712^2] [^2+^352] 

6 r ^ , n> l2 



7714 — [Ai + 77iii?i] [Ai + m/^B 



Here ttii, 7774 are the magnetizations of the two different representative sites in the core and 7772, 7773 
are the magnetizations of the two different representative sites in the shell. The coefficients in Eq. 
® are 

An = cosh ( J„V) F {x) \x=o 
Bn - sinh(J„V)F( X ) \x=0 

where V is the usual differential operator in the differential operator technique and n — 1,2,3. 
The function is defined by 

F{x)^ JdH,P{H,)f{H,,x) (5) 

where 

f {Hi, x) =teiiih{/3x + /3Hi) . (6) 

In Eq. /3 — l/{kBT) where kg is Boltzmann constant and T is the temperature. The effect of 
the exponential differential operator to a function is given by 

exp (aV)F {x) ^ F {x + a) (7) 

with any constant a. DA will give the results of the Zernike approximation [47' for this system. 
With the help of the Binomial expansion, Eq. ([3]) can be written in the form 

4 12 1 

mi = ^'^'^'^Ki{i,j,k,l)m\mlm^Tn[ 



i=0j=0k=0l=0 
12 2 



7772 = X!X!X!^2 {i, j, k) m\ 



7 h 
7772 7773 



i=0j=0k=Q 



2 2 2 

7773 = {i, j, k) m\m'2^^ 

i=0j=0k=0 
6 2 

7774 = if 4 (7, Z) 7775^7774 

1=0 1=0 



(8) 



where 



Ki {i,j,k,l) 



4 \ / 1 \ / 2 

* / \ / \ ^ 

1 \ / 2 \ / 2 

2 W 2 \ / 2 



1 ^ ^5-^-i^3-j-fc^j+;^j+/c 



K2{i,j,k) -(:)(:)( : i ^-^-'^Ai-si+^s; 



K3{^,J,k) = { ; { j [1 ]Ar-^Ai-Wr^Bl 
Ki{i,l) = (!)(?) Al'"''B\+' 



2 ^3 ^2 -"3 

(9) 
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These coefficients can be calculated from the definitions given in Eqs. ^ and ([7]). 

For a given Hamiltonian and field distribution parameters, by determining the coefficients from 
Eq. ([9]) we can obtain a system of coupled non linear equations from Eq. ([8]) , and by solving this 
system we can get the magnetizations mi,i — i,2,3,4. The magnetization of the core {rric) and 
shell (ttIs) of nanowire, as well as the total magnetization (tot) can be calculated via 



TOc = - (6mi + TO4) , TOs = — (6to2 + 6TO3) , 
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(6toi + 6to2 + 6TO3 + TO4) (10) 



Since in the vicinity of the critical point all magnetizations are close to zero, we can obtain 
another coupled equation system for determining this critical point by linearizing the equation 
system given in Eq. ([8]), i.e. 

A.m^Q (11) 



where 



( Xi(l, 0,0,0) 

^2(1,0,0) 
^3(1,0,0) 

if4(l,0) 



i^i(0, 1,0,0) JCi(0,0,l,0) ifi(0,0,0,l) \ 



V 



^2(0,1, 
^3(0,1, 





0) 
0) 



^2(0,0,1) 

^3(0, 0,1) 








if4(0, 



} 



TO 



I mi \^ 

TO2 

TO3 

\ 1714 J 



(12) 



(13) 



As discussed in [33^, the matrix A 
- J3 then we can conclude that the 



Critical temperature can be determined from det{A) — 0. 
given in Eq. (|12p is invariant under the transformation J3 — >■ 
system with ferromagnetic (J3 > 0) core-shell interaction has the same critical temperature as 
that of the system with anti- ferromagnetic (J3 < 0) core-shell interaction (with the same | J3I) for 
certain Hamiltonian and magnetic field distribution parameters. Although this discussion has been 
made for the system with zero magnetic field in [33 , this conclusion is also valid for this system, 
because of the symmetry of the magnetic field distribution. The equation det(A) = is invariant 
under the transformation J3 — > — J3 for the nanowire with trimodal magnetic field distribution. 



4 Results and Discussion 

In this work we discuss the effect of the random magnetic field distribution on the phase diagrams 
of the system. Since the phase diagrams are the same for the J3 > and J3 < with the same 
IJ3I, we focus ourselves on the case J3 > 0, i.e ferromagnetic core-shell interaction. We use the 
scaled interactions as 

Ji = J, '^n^y, n = 2,3. (14) 
Let us start with bimodal magnetic field distribution. 



4.1 Bimodal Distribution 

This distribution is given by Eq. ^ with p = 0. It is a well known fact that, in this distribution, 
increasing Hq values drag the system to the disordered phase, due to the random distribution of 
the magnetic fields zLHq on the lattice sites. On the other hand, the interactions J„, (n = 1,2,3) 
enforce the system to stay in the ordered phase. Thus a competition takes place between these 
two factors, in addition to the thermal agitations. 

The phase diagrams in [kBTc/ J, Hq/ J) plane can be seen in Fig. [2] for different r2 and 
values. We can see from Fig. [2] that the general effect of increasing Hq is to decrease the critical 
temperature, as expected. Besides, another expected result is that the higher r2 or r3 values 
has the higher Tc values at a certain Hq value which can be seen in Fig. [2] But we can see 
in Fig. [Ha) that, the latter observation does not hold for r3 = case. In the absence of the 
interface interaction between the core and shell, the phase diagrams in (kBTc/ J, Hq/ J) plane are 
the same for all r2 < values. In order to clarify this situation, we plot the curves of the core 
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magnetization, shell magnetization and the total magnetization of the system in Fig. [31 As seen 
in Fig. [21 shell magnetization falls to zero at a temperature lower than that of the core for the 
r2 — 0.8, 1.0. For these values, the core determines the critical temperature of the system. But for 
the value of r2 — 2.0, the critical temperature of the shell is higher than the critical temperature 
of the core. Then we can say that shell of the nanowire has a lower critical temperature than 
that of the core, up to a certain r2 — value, after that value, the critical temperature of the 
shell becomes higher than the critical temperature of the core i.e. the shell begins to determine 
the critical temperature of the system. This means that for the values that satisfies r2 < rl^, 
the phase diagram of the system does not affected by the changing r2- Although the r2 — 1.5 
curve in Fig. [H^a) is different from other r2 curves corresponding to r2 < 1.5 for the higher Hq/J 
(which are close to the Hq/J = 3.0) values, we can say that this r2 value satisfies r2 > 1.0. 
This fact can be explained by taking a look at the geometry of the system. In the absence of 
the interaction between the core and shell, both of the two distinct shell spins have four nearest 
neighbors (which have magnetizations labeled as 7712,7713) and the two distinct core spins have five 
and eight nearest neighbors, respectively (which have magnetizations labeled as 7711,7(74). Thus we 
can say that for equal exchange interactions of the core and shell (72 = 1), the core has higher 
critical temperature than that of the shell. Consequently, 72 value that providing the equality of 
the critical temperatures of the core and shell has to be 72 > 1.0. In summary, we can conclude 
that, for the ^3 = phase diagrams in (ksTc/ J, Hq/ J) plane, all diagrams are the same up to a 
certain 73 > 1.0 value, and when 72 = 7'2, the diagrams begin to separate from the diagrams which 
are for the r2 < 7*2 in the higher Hq/J region. Furthermore, as r2 grows then the phase diagrams 
become completely separate from the diagrams for the 72 < 7'2. 

We can see from Fig. [51 that, increasing 72 and 73 values make the ferromagnetic region wider 
in (kBTc/ J, Hq/ J) plane, since increasing 72 and 73 values rises the absolute value of the lattice 
energy coming from the spin-spin interaction, which must be overcomed by both thermal agitations 
and random magnetic fields for the phase transition from an ordered phase to a disordered one. 




1 2 3 4 5 01234567 012345678 



Hq/J Hq/J Hq/J 



Figure 2: The phase diagrams of the nanowire with bimodal random field distribution in the 
{ksTc/ J, Hq/ J) plane for different 72,73 values. 

Now the question is: what is the magnetization behavior for that different 72, 7*3 and Hq values 
for bimodal magnetic field distribution. Since the magnetization behavior with temperature is well 
known for the nanowire without any random magnetic field and due to the number of parameters 
in this problem and also the restrictions on the length of the paper, we restrict ourselves to 
this question: what is the ground state magnetization behaviors with Hq form different 7*2,73 
values? We can see the answer of the question from Fig. [H which is the variation of the core 
ground state magnetization(77i^), shell ground state magnetization (mf ) and the total ground state 
magnetization (m^) (where the superscript g stands for the ground state) with Hq. This figure 
corresponds to the phase diagrams given in Fig. [21and the temperature is chosen as fc^T/ J = 0.001 
which can be regarded as ground state for the system. As seen in Fig. [IJg), in the absence of 
the core-shell interaction (i.e. 73 = 0), two more ground state magnetization appears apart from 
the magnetizations 777^, = 0.0 and 777^ = 1.0 which are correspond to the disordered state and 
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Figure 3: Variation of the magnetization with temperature for some selected r2 parameter value of 
the nanowire with bimodal random field distribution. Fixed parameter values are rs = 0.0, Hq/ J = 
1.5. 

ordered state, respectively. These two more ground state magnetization values are m|, = 0.3684 
and — 0.6316 and can be regarded as two different partially ordered phases of the nanowire 
under the bimodal magnetic field distribution. Since the core has the ground state magnetization 

= 1.0 regardless of the value of r2 up to the value of Hq/J = 3.0 ( Fig. HJa)), that partially 
ordered phases seen in Fig. IH^g) has to come from the mismatch in magnetization behaviors with 
Hq in the core and shell. It can be seen from Eq. ([T0| that the core contributes to the total 
magnetization with the ratio 7/19 while this ratio for the shell is 12/19. Because of the difference 
of the regions which are core and shell in ordered and disordered phases with Hq , this two partially 
ordered phases appears in the regions seen in Fig. IDJg). For example, for the = 0.0, r2 = 2.0, 
the core has ground state magnetization = 1.0 up to Hq/J = 3.0 and the shell has ground state 
magnetization mf = 1.0 up to Hq/J — 4.0 then a partially ordered phase = 0.6316 = 12/19 
appears in the intersection of this two region which is 3.0 < Hq/J < 4.0. 

The presence of the interaction between the core and shell makes this situation slightly different. 
When ra ^ 0.0 then the system becomes whole instead of separate core and shell as in the case 
ra = 0.0. Then the core and the shell can develop partially ordered phases separately, thus new 
partially ordered phases can appear in the system. As seen in Fig. HJe), shell can develop two new 
partially ordered phases (other than phases that corresponds to ordered mf = 1.0 and disordered 
mf = 0.0 phases) which have magnetizations mf = 0.5 (r2 = 0) and mf = 0.9835 (r2 — 0.5, 1.0) 
due to the random field and presence of the interaction with the core. This makes the ground state 
magnetization of the system as in Fig. lU^h) with the same reasoning explained just above, with 
two partially ordered phases namely = 0.6842 (r2 = 0) and m|, = 0.9896 (r2 = 0.5, 1.0). The 
increasing interaction strength between the core and shell increases the magnetization of that two 
partially ordered phases that has magnetization near to 1.0 of the shell and then the ground state 
magnetizations of the system become as seen in Fig. IH^i). 

More detailed investigation on that partially ordered phases can show that, indeed a few other 
partially ordered phases exist in the shell even in the core. This fact can be seen in Fig. \E\ As in 
Fig. lUJa) for = 1.0, with little 7-2 (e.g. r2 = 0.3), rising Hq can induce a partially ordered phase 
in the core which disappears with rising r2 (e.g. r2 = 0.5). The same thing also be valid for the 
shell as seen in Fig. [Sfb) and thus whole system which can be seen in Fig. [S](c). 

4.2 Trimodal Distribution 

Let us investigate the phase diagrams in the (ksTc/ J, Hq/ J) plane for the distribution which is 
given by Eq. ([2]) for the p ^ 0. We can see from the random magnetic field distribution that, 
when p approaches to p = 1 from the value p = 0, the phase diagrams given in Fig. [2] have to 
evolve to the phase diagrams of the system with zero magnetic field, i.e. parallel lines to the 
Hq/J axes in that plane. This evolution can be seen in Fig. [B] We can see from the Fig. IHlthat 
(e.g. r2 = 1.0, = 1.0 curve in order of Fig. [6] (a),(c),(e),(g)) , rising p from the p ~ Q to the 
p = 1 first opens phase diagrams on the right hand sides of them, then creates a tail on that side 
which extends to the high magnetic fields then rises that tail to the level of the pure system's 
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critical temperature. The same effect can be also performed by the rising r2 (or r^) in the trimodal 
distribution for little p values, since that effect results from the competition between the spin-spin 
interaction which tends system to ordered phase and randomly distributed magnetic field. This 
can be seen for instance in Fig. |6l[c). We also note that, during this evolution process a phenomena 
called reentrant can be seen in the diagrams. 

The details of the evolution process of the phase diagram with p in the {ksTc/ J, Hq/ J) plane 
for fixed r2 — 1-0, — 1.0 values can be seen in Fig. [T] We can see from Fig. [7] that, at zero 
temperature the transitions from the ordered phase to a disordered one occurs at integer Hq/J 
values. 

The ground state magnetizations appear more complicated forms for this distribution than 
the bimodal distribution. This can be seen from Fig. |S] which are corresponding ground state 
magnetizations to the phase diagrams given in Fig. |n^a),(c),(e),(g). Construction of the total 
ground state magnetization from the ground state magnetizations of the core and shell is the same 
as explained in the bimodal distribution section, since trimodal distribution does not change the 
contribution ratios of the core and shell to the total magnetization. But we can see from Fig. |8] 
that, trimodal distribution induce higher number of partially ordered phases to the system both 
core and shell than the bimodal distribution. This causes higher number of partially ordered phases 
of the whole system. It is not meaningful to give all the ground state magnetization values of that 
ordered phases seen in Fig. [51 since that values depends on the parameters r2,ry,,p and Hq. But 
we have to say that rising r2 for fixed = 1.0, shifts the partially ordered phases to the right of 
the Tot ~ Hq/J plane, means that increasing ^2 destroys partially ordered phases with beginning 
in low Hq/J. We can also see from Fig. |S]that the same effect can be done by rising p, i.e. also 
rising p destroys the partially ordered phases beginning from the low Hq/J. Thus we can say that 
rising r2 for other fixed parameters or rising p for other fixed parameters induces same changes 
on both phase diagrams and the ground state magnetizations. This similarity is expected because 
rising r2 means that interactions between spins in the shell getting stronger, then after some value 
of 7^2 this interactions can overcome the effect of the random magnetic field distribution. On the 
other hand rising p means that the falling strength of the randomness, then system can escape 
from that partially ordered phases which are formed by randomness effects. 

5 Conclusion 

The effect of the random magnetic field distribution on the phase diagrams and ground state 
magnetizations of the Ising nanowire investigated. The evolution of the phase diagrams of the 
system with bimodal magnetic field distribution to the pure one (i.e. system with no magnetic 
field) investigated as the distribution parameter of the magnetic field goes from p = to p = 1. In 
the bimodal distribution as a special case of the trimodal distribution with p — 0, some partially 
ordered phases seen for the system as a result of the inspection of the ground state magnetizations. 
This complicated situation in the ground state magnetizations comes from the competition between 
the randomness of the magnetic field which drag the system to the disordered phase and the 
interaction strengths between the spins which try to hold the system in an ordered phase. 
Our findings about the phase diagrams in (fcsTc/J — Hq/J) plane as follows: 

• All closed phase diagrams intersect Hq/J axes at integer Hq/J values. This means that at 
zero temperature, system has ferromagnetic ordering up to the some value of Hq which is 
multiple of the core interaction strength J. 

• In the absence of the core-shell interaction {rj, = 0.0) the phase diagrams for different r2 
values are the same up to some r2 < value where r2 > 1.0 for the bimodal distribution. 
This fact comes from the difference between the total interaction energy of the core and shell. 
When r2 reaches the value rl^ , shell starts to determine the critical temperature of the system 
and thus phase diagrams for r2 > ^2 become separate. 

• Phase diagrams of the system with trimodal magnetic field distribution getting open on the 
high Hq/J region when p is sufficiently large. This is due to the decreasing randomness. 

Lastly our findings about the ground state magnetization values as follows: 
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• System with bimodal magnetic field distribution have a few partially ordered phases. 

• This number of partially ordered phases getting higher in the trimodal distribution. 

• Rising p and also rising r2 or destroys this partially ordered phases. 

These partially ordered phases requires more detailed investigation in order to obtain more 
precise relations of the numerical values of that ground state magnetizations and the length of a 
certain ground state magnetization in the (m — Hq/J) plane, with random magnetic field distribu- 
tion parameters. At the same time these partially ordered phases means that, the phase diagrams 
in {ksTc/j — Hq/J) plane have some first order transition lines. But since one can not calculate 
the free energy of the system within this formulation, we can not explicitly determine these first 
order transition lines. 

Altough these questions waiting for being answered, another important question is the eff'ect 

of the continuous random magnetic field distributions (e.g. Gaussian distribution) to the phase 
diagrams of the Ising nanowire. This will be our next work. We hope that the results obtained in 
this work may be beneficial form both theoretical and experimental point of view. 
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Figure 4: Variation of the ground state magnetizations with Hq/J which corresponds to the 
phase diagrams seen in Fig. [2] for the system with bimodal magnetic field distribution. Each 
magnetization calculated at a temperature ksT/J = 0.001 which can be regarded as ground state. 
The letters which are close to the related curve represent the different r2 values as (A) r2 — 0.0, 
(B) r2 = 0.5, (C) r2 - 1.0, (D) = 1.5, (E) ra = 2.0 
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Figure 5: Variation of the ground state magnetizations for selected values of r2 for the system with 
bimodal magnetic field distribution. Each magnetization calculated at a temperature ksT / J = 
0.001 which can be regarded as ground state and for the parameter value rj, — 1.0. 
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Figure 6: The phase diagrams of the nanowire with trimodal random field distribution in the 
{ksTc/ J, Hq/ J) plane for different r2,r^,p values 
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Figure 7: The phase diagrams of the nanowire with trimodal random field distribution in the 
{kBTc/ J, Hq/ J) plane for different p values. Fixed values of the diagrams are r2 = 1.0 and rs = 1.0. 
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Figure 8: Variation of the ground state magnetizations for selected values of r2 and p for the 
system with trimodal magnetic field distribution. Each magnetization calculated at a temperature 
ksT/J = 0.001 which can be regarded as ground state and for the parameter value = 1.0. The 
letters which are close to the related curve represent the different p values as (A) p = 0.0, (B) 
p = 0.3, (C) p = 0.5, (D) p = 0.7, (E) p = 1.0 
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